Unconventional superconducting phases on a two-dimensional extended Hubbard 

model 
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We study the phase diagram of the extended Hubbard model on a two-dimensional square lattice, 
including on-site (f/) and nearest-neighbor [V) interactions, at weak couplings. We show that the 
charge-density- wave phase that is known to occur at half-filling when 'iV > U gives way to a dxy-w&ve 
superconducting instability away from half-filling, when the Fermi surface is not perfectly nested, 
and for sufficiently large repulsive V and a range of on-site repulsive interaction U. In addition, when 
nesting is further suppressed and in presence of a nearest-neighbor attraction, a triplet time-reversal 
breaking (px + ipy)-wave pairing instability emerges, competing with the dx2^y2 pairing state that 
is known to dominate at fillings just slightly away from half. At even smaller fillings, where the 
Fermi surface no longer presents any nesting, the {px + ipy)-wa.ve superconducting phase dominates 
in the whole regime of on-site repulsions and nearest-neighbor attractions, while d^y-pairing occurs 
in the presence of on-site attraction. Our results suggest that zero-energy Majorana fermions can 
be realized on a square lattice in the presence of a magnetic field. For a system of cold fermionic 
atoms on a two-dimensional square optical lattice, both an on-site repulsion and a nearest-neighbor 
attraction would be required, in addition to rotation of the system to create vortices. We discuss 
possible ways of experimentally engineering the required interaction terms in a cold atom system. 
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I. INTRODUCTION 

An extended Hubbard model is generally employed as 
a theoretical framework of screened electronic interac- 
tions and regarded as a prototypical scenario for rich 
quantum phases in condensed matter physics^"— . In a 
onc-dimcnsional chain, for instance, an extended Hub- 
bard model, including an on-site and nearest-neighbor 
interactions, presents correlated phases associated with 
the ratio of the two interactions^. Recent indentifi- 
cation of the bond-charge-density-wave instability be- 
tween charge and spin density- wave phases, at weak 
and strong interactions, completes the phase diagram 
of this model^"— . Meanwhile, in a two-leg ladder, a 
checkerboard charge-ordered state has been proposed for 
all fillings between quarter and half, with on-site and 
nearest-neighbor repulsioi:^^^. Its application to the cou- 
pled quarter-filled ladders with coupling to the lattice has 
recently been studiedi^ to explain the spin gaps in the 
NaVzOs materialii^. 

On a two-dimensional lattice, the extended Hub- 
bard model has been considered a paradigmatic model 
to search for possible unconventional superconducting 
phases since the discovery of high-temperature supercon- 
ductivity in the cupratea^^. Although a nearest-neighbor 
repulsion between electrons suppresses non-s-wave pair- 
ings tcndenciesii, it is generally believed that the nesting 
of the Fermi surface plays a key role in driving uncon- 
ventional pairing under purely repulsive interactions at 
weak couplingai^. In the proximity of density-wave or- 
der, for instance, a chiral d-wave state has been found for 
an extended Hubbard model on both triangular and hon- 
eycomb lattice&i^r— . Furthermore, following the recent 



experimental realization of a two-dimensional Kagome 
lattice for ultracold atoms2^, the phase diagram of the 
extended Hubbard model on a Kagome lattice has been 
established in the vicinity of van Hove fillings^^. It was 
shown that a possible p-wave charge and spin bond order 
can be triggered in the presence of a nearest-neighbor re- 
pulsion for van Hove fillings, then giving way to a /-wave 
superconducting phase when slightly doped awaj^. 

It has recently been proposed that under a nearest- 
neighbor attraction and on-site repulsion a singlet (p + 
ip)-wave pairing emerges on a honeycomb latticG^^. 
Along with the result that Majorana fermions can be 
generated as a zero-energy mode in the excitation spec- 
trum of a half-quantum vortex in a {p + ipywave 
superconductor—, this indicates the possibility of creat- 
ing Majorana fermions in graphene in the presence of 
a magnetic field. However, a functional rcnormalization 
group study has shown that for a honeycomb structure, 
the /-wave pairing is preferred and is stabilized by intro- 
ducing a next-nearest-neighbor attractior^i^. Although 
triplet p-wave superconductivity has been proposed for 
an extended Hubbard model on a square lattice for purely 
repulsive interactions^i, the physics of long-range interac- 
tions, with possibly terms leading to competing instabili- 
ties, has been studied recently, leading that the existence 
of p-wave pairing is still an open questio n^^'^^ . 

In this paper, we study the phase diagram of an ex- 
tended Hubbard model via a functional rcnormaliza- 
tion group (fRG) approach^Si^j including on-site U and 
nearest-neighbor V interactions, on a two-dimensional 
square lattice. The total Hamiltonian can be written as, 
H = Hq + i/int , with the noninteraciting and interacting 
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FIG. 1: (Color online) Fermi surfaces of a square lattice with 
nearest-neighbor hopping, at different chemical potentials fj, 
(in units of t). Below, the phase diagram versus /i for U — 1 
and V = 0. The Fermi surface patches used in this study are 
illustrated in the Fermi surface at /i = 0. 
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FIG. 2: (Color online) The phase diagram (a) at half-filling 
fi = 0, (b) at /i = —0.095, parameterized by on-site U and 
nearest-neighbor interaction V. The interaction terms U and 
V are in units of t throughout this paper. 



parts, 

Ho = -tJ2 (cLcja +H.C.) -fij^n,, (1) 

{ij),a I 

ni-fU^l + V" ^ TijTij , (2) 

respectively, where (ij) represents nearest-neighbor pairs 
of sites, rii = ni-f + riii = J2a '^ia^ia ^ is the chemical 
potential. Without the nearest-neighbor interaction, the 
phase diagram in a fRG analysis is well developed in the 
limit of weak couplings^. For on-site attraction, s-wave 
superconductivity (s-SC) dominates at all fillings except 
for half-filling. The phase diagram for on-site repulsion, 
U ~ 1^ versus /i is sketched in Fig. [TJ In the vicinity 
of half-filling, the spin-density wave (SDW) dominates 
due to the strong nesting of the Fermi surface. With 
slight doping, c?2.2^,y2-SC emerges in a small regime of 
/i just away from half-filling. When the magnitude of 
/i is further increased, no instability develops up to the 
point where wc stop the RG flows, at energy cutoffs lower 
than 10~^t. Kohn-Luttinger (KL) instability^, with ex- 
tremely low critical temperatures, is expected to occur in 
this regime. 

In presence of a nearest-neighbor interaction V , the 
phase diagram is much richer. We find that dxy- and 
{Px + ipy)-"wave pairing superconducting states develop. 
In the proximity of the charge-density-wave (CDW) or- 
der, dajy-wave pairing emerges from the CDW instability 
with U > and V > 0. When the nesting of the Fermi 
surface is decreased, a time-reversal symmetry breaking 
(Px + ipi/)-SC arises from the (ij.2_|_,y2-SC with a suffi- 
cient large nearest-neighbor attraction. When nesting 
is completely suppressed, (px + ipy)-SC dominates in the 
whole regime of U > and V < 0, and the d^y-pairing 
is only triggered with the help of an on-site attraction. 
Using a symmetry argument, we show that appearance 
of {px + ipy)-SC on a square lattice for < is generic 



and robust due to the underlaying lattice structure, and 
can be used to create a zero mode Majorana fermion in 
the presence of a magnetic field. 

This paper is organized as follows. In Sec. II, we study 
the phase diagrams at and near half-filling (/i = and 
—0.095). We will show that our fRG results are consistent 
with the previous studies at the half-filling. In Sec. Ill, 
we show the phase diagrams for chemical potential /i = 
—0.5, —1 and —2 and study the possible unconventional 
SC states driven by the nearest-neighbor interaction. We 
discuss our results and conclusions in Sec. IV. 



II. AT AND NEAR HALF FILLING 

Starting with the bare Hamiltonian, Eq. ([2]), we fol- 
low standard fRG procedure integrating out high-energy 
modes, decreasing the energy cut-off A. The four-fermion 
terms in the resulting effective Hamiltonian are written 

in the form g(fci, fc2, fca, A)V'l(fci)'0]j(fc2)V'/3(fc3)'0Q(fci + 
^2 — fea), in momentum space, with k = {kx,ky) and 
spin indices a, (3. In previous studies, the RG equa- 
tions for models with SU{2) and C/(l) symmetries have 
been systematically studied by Fermi surface discrete 
patch-approximatior*^. In this paper, to preserve the 
particlc-holc symmetry of the non-interacting Hamilto- 
nian at half-filling, all phase diagrams are obtained by the 
configuration of the Fermi surface patches illustrated in 
Fig.[TJ By integrating out high energy degrees of freedom 
and neglecting self-energy corrections, the RG flows of all 
couplings versus the decreasing running energy cutoff A 
are computed. Before the system flows into the strong 
couplings regime, we truncate the RG process when the 
absolute magnitude of one of the couplings reaches ~ 30i. 

To determine the dominant instability, we decompose 
specific four-fermion interaction terms in the Hamilto- 
nian as X^ife p ^op(fc,P, A)C'^C'p, with Ok a bi-fermion 
operator for the order parameter (op) of SC, CDW, SDW 
or Pomeranchuk instability. Then, for a given order pa- 
rameter channel, wc further decompose, Vop(fe,p, A) = 
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FIG. 3: (Color online) Phase diagrams for (a) n = —0.5, (b) /j. — —1 and (c) /j. — —2, parameterized by U and V. The 
transparent regime means no instability is found before we stop the RG process, at which point the energy cutoff is lower than 
W^^'t. Form factors obtained from decoupling of pairing channels into bi-fermions during the fRG flows are illustrated for (d) 
dxy-wave (e) (px + ipy)-wsLve pairings. 



X;j<p(A)/*;(A;, A)/*p(p, A), in normal modes, with i 
a symmetry decomposition index. The leading insta- 
bility can be determined by the most minimum eigen- 
value w™™(A) (largest magnitude), and the correspond- 
ing symmetry of the instability is given by the form factor 

/-i"(fc)33^. 

At half-filling, the phase diagram has been studied ex- 
tensively by several methoda^ '^^i'^^ . In our fRG analy- 
sis, we include the regime of negative interactions and 
obtain the phase diagram in Fig. [2^, parameterized by 
on-site interaction U and nearest-neighbor interaction V 
(in unit of t). The phase boundary between SDW and 
CDW is at J7 ~ iV for U,V > 0, consistent with known 
results from previous studies. For the line of I^ = and 
U < 0, we find that CDW and s-SC are degenerate, also 
in agreement with results in the literature^. However, 
the degeneracy is delicate and broken by introduction of 
a nearest-neighbor interaction: a slight nearest-neighbor 
attraction drives the system to a s-SC instability, instead 
of a CDW instability. 

We also find that for a sufficiently large nearest- 



Ill. DOPED SYSTEMS 



Here, we increase doping, decreasing nesting of the un- 
til the density-wave instability no longer occurs. Then, 
as shown as Fig. [3^, the dxy-wave SC arises from the 
CDW instability in the regime oi V > and U > 0. 
The form factors of the c?^j,-SC in our fRG analysis is 
plotted in Fig. [3ji. Although the dxySG phase has been 
proposed for the purely repulsive modela^I^iSS^ the dxy- 
SC instability we find develops only with an appropriate 
nearest-neighbor repulsion and also tied with the prox- 
imity to nesting of the Fermi surface. It is only when 
nesting is barely suppressed so that CDW is no longer 
dominant, but CDW fiuctuations are still expected to be 
strong, that this d^y-SC phase emerges. 



neighbor attraction, a d^2_|_j,2-SC is triggered, even dom- 
inant over the s-SC in the regime of [/ < 0. This d-wave 
pairing is linked with the nesting of the Fermi surface. In 
other words, if the nesting effect is suppressed, the s-SC 
is the dominant instability for a generic Fermi surface 
in the regime of /7 < 0. Furthermore, without nesting 
effects, a p-wave, with lower angular momenta than d- 
wave, will eventually dominate in the regime of < 
and U > 0. However, negative V combined with nesting, 
leads to dx2+y2-SC. 

By slightly doping away the half-filling, the dx2-^-y2-SG 
overcomes the spin-density-wave instability that domi- 
nates for on-site repulsion. The phase diagram parame- 
terized by U and V is illustrated in Fig. Wp- The SDW 
is suppressed to the small regime between CDW and 
d.j.2^y2-SC in the phase diagram. The degenerate line, 
U < and V = 0, mentioned above, is dominated by 
s-SC instability since the Fermi surface is not perfectly 
nested. However, the CDW is still dominant in the over- 
all regime otV > 0. 

I 

As compared with the phase diagram at /.t = —0.5 
in Fig. [5^, the dxy-SC instability moves to the region 
of larger positive V or negative U when doping is in- 
creased, for fi = —1, as shown in Fig. |3}d. By looking 
at the RG fiows of different couplings, we notice that 
the different behavior for the two fillings, /i = —0.5 
and /i = — 1, are mainly coming from the effect of the 
nesting vertices, that is, couplings (?(fci, ^2, fca, fc4) with 
fci + fc2 - fcs - A;4 = (0,27r) or (27r,0) or {2tt,2tv). For 
/i = — 1, there are fewer nesting vertices flowing into large 
values as compared with the /i = —0.5 case. Without 
the benefit of nesting vertices, the instability of d^y- and 
dx2j^y2-SC are suppressed and eventually do not develop 
by the time we stop the RG fiow. As a consequence, when 
the nesting is completely destroyed, in heavily doped 
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cases, there is no instability found for purely repulsive 
interactions {U > and V > 0), a.s shown in Fig. 

The triplet p-SC phase tells an opposite story. In the 
regime of < 0, under the influence of Fermi surface 
nesting, dx2_^_y2-SC is dominant, as shown as Fig. [3^. 
However, decreasing nesting suppresses the d-wave pair- 
ing and the p-SC emerges for large nearest-neighbor at- 
traction, as shown in Fig. ^p. When the nesting is no 
longer present at all, p-SC dominates the regime of F < 
and U > 0, as illustrated as Fig. [3J;. The emergence of 
a triplet p-SC has a clear physical picture in a square 
lattice. In the presence of a nearest-neighbor attraction 
V < 0, electronic pairings are triggered. To lower energy, 
a s-wave pairing with the lowest angular momenta would 
be favored. However, an on-site repulsion suppresses the 
s-wave pairing so that p-wave, with the second lowest 
angular momenta, is preferred. 

Furthermore, the four-fold symmetry on a square lat- 
tice also indicates that the Px- and py-wave instability 
channels must be degenerate. The degeneracy of two p- 
wave superconducting instabilities is indeed found in our 
fRG results and the associated form factors are plotted 
in Fig. [3^. Since the p^- and p^-wave superconducting 
states are degenerate, the gap functions will be given 
as a linear combination of the two order parameters. By 
constructing the gap function as = Ap^ (fc) -|-wAp (fc) 
with the complex coefficient v containing a possible rela- 
tive phase, the condensation energy in the standard BCS 
equation is obtained by the difference in energy between 
the superconducting and normal states, and is given by 
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V4 + I^fcl 



(3) 



where is the dispersion relation of the non-interacting 
Hamiltonian. The second term of Eq. (|3]) is maxi- 
mized when V is purely imaginary, hence the time-reversal 
breaking pairing symmetry px -t- ipy is the energetically 
favored one^Si^i^i^. Physically this is reasonable, since 
this choice of the phase guarantees that a gap forms ev- 
erywhere along the Fermi surface, lowering the ground- 
state energy. 



wave superconducting state has been proposed in these 
systems^. However, the mechanism discussed here puts 
some constraint for the emergence of the Px + ipy-SC 
phase: in order to develop the px -\- ipy-SC, long-range 
attraction and on-site repulsion is needed, as well as low 
density of fermions to avoid nesting of the Fermi sur- 
face, while large enough densities such that the Fermi 
energy is a large scale (t > U, V) to justify the validity of 
fRG results. Another possible way to manifest a medi- 
ated attraction is through another species of fermions in 
a Fermi-Fermi atom mixture^. By introducing an inter- 
species interaction in Fermi-Fermi mixtures, an effective 
interaction for one species of fermions can be obtained 
by tracing out the other species. This can be justified, 
for example, if the one of the species has a much smaller 
effective mass than the other. In this caso^, the me- 
diated long-range interaction is found to decay rather 
rapidly and can be approximated by an effective on-site 
and nearest-neighbor interactions. By tracing out one 
species with low electronic density, an effective nearest- 
neighbor attraction can be obtained. Together with a 
bare hard-core on-site repulsion, it may provide the re- 
quired conditions for the creation of the time-reversal 
breaking px -\- ipy-w&ve pairing. 

In conclusion, we study the phase diagram of an 
extended Hubbard model, including a on-site U and 
nearest-neighbor V interactions, on a two-dimensional 
square lattice. In the proximity of charge-density-wave 
order, the dxy-SC overcomes the CDW, dominating in 
the regime of J7 > and V > 0, from our fRG analysis. 
Accompanying the destruction of Fermi surface nesting, 
a time-reversal breaking {px -\- ipy)-wa,ve superconduct- 
ing state arises in the regime of V < 0. Our results 
indicates that, without nesting, the {px 4- ipy)-SC on a 
square lattice under a nearest-neighbor attraction is the 
generic behavior due to the underlaying lattice structure, 
and can be used to create a zero mode Majorana fermion 
in the presence of a magnetic field. 
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IV. DISCUSSION AND CONCLUSION 

In atomic Bose-Fermi mixtures, an effective attraction 
between fermions can be mediated by fluctuations of the 
Bose-Einstein condensate of the bosons^"—. In the pres- 
ence of the mediated long-range attraction, the Px + ipy 
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